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1.0 PROGRAM DESCRIPTION 


The program BUCLAP2 uses minimum energy principles to do an elastic stability analysis of a 
curved or flat, laminated, rectangular long plate subjected to combined inplane normal and shear 
loads. 

The structure of the program can be understood from the fundamental problem which is being 
solved. The fundamental problem is: 


|S„| -0 (1) 

where S a is the stiffness matrix of the plate under consideration with the appropriate boundary 
conditions applied for a given load, a. This problem is classified as a general nonlinear eigenvalue 
problem. The elements of S a are transcendental functions of the external loading, a, and the axial 
half-wave length qf buckling (a/m). 

For any chosen number of axial half-waves (a/m), the lowest level of external loading at which 
equation 1 is satisfied is the buckling load of the structure. This load is determined by an iteration 
procedure. A series of (a/m) values are investigated and the lowest of all buckling loads then gives 
the critical load of the structure 

The above-mentioned iteration procedure may be described as follows: For a given (a/m) value 
there is an associated interval, og < a < a u , of interest (this interval is either calculated or specified 
by input). A modified bisection search technique is applied to the interval. So for the first iteration, 
a equals og + f(a u - a g), (0 < f < 1). For this load, a, the stiffness matrix for the plate is calculated 
using equation A.53 of reference 1 . The elemental matrix is merged into S a by applying the 
boundary conditions and stiffnesses of the edges of the plate. 

The resulting S G matrix is complex Hermitian. The choice between subinterval [ag, a] and 
[a, n u ] that is of further interest is dictated by whether the determinant is positive definite, or 
not, at a. If positive, the right subinterval is chosen as the new interval; if not, the left. The 
objective is to reduce the length of the interval being investigated to a suitable tolerance. As the 
quantity a (which is being sought) is changed, S ff changes in a complex manner. This character- 
istic affects the structure of the program profoundly. 



1.1 PROGRAM FUNCTION 


The program consists of three primary overlays with one primary having three secondaries. 

Overlay Function 

0.0 Monitors basic control 

1 .0 Inputs and processes data 

2.0 Monitors control through load calculation 

2. 1 Calculates wave numbers for minimum search 

2.2 Calculates upper bound loads for list of wave numbers 

2.3 Obtains minimum load of plate for list of wave numbers 

4.0 Summarizes load calculation at minimum value 

The process of producing the plate merged stiffness matrix S a is done in the routines DATIN, 
DC, MERGEC and STOREC. A discussion of this process follows. The plates may be visualized 
as consisting of two nodes, each node with four degrees of freedom, w, 0, v, u: 

w - out-of-plane displacement 

0 - rotation about plate sides 

v - inplane displacement in lateral direction 

u - inplane displacement in longitudinal direction 

- stiffness matrix of half-plate 

S - stiffness matrix of “full plate” 

The merged stiffness matrix for a plate P is produced by generating the stiffness matrix for a 
half-plate with width one-half of P. This matrix is used in producing both the left part and right 
part of P. The left part consists of the matrix for half-plate P^ and the boundary conditions and 
nodal stiffnesses at the left edge of P. The right part consists of the matrix for half-plate and the 
boundary conditions and nodal stiffnesses of the right edge of P. 
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Discussion of Merge procedure follows: 


Given Xj, Xf, then • X^ = Xf or = Xf • Xj"* 


£ r 

1 2 3 4 1 2 3 4 



M R 

Merged Matrix Before Row Deletion 


Both matrices (dashed and dotted) have the same members before deletion. The deletion of a 
freedom for a plate side is accomplished by deleting the row and column of the merged matrix 
corresponding to the side and freedom (i.e., if w = 0 for side L then column corresponding to Lj is 
deleted and row Lj is also deleted). 
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Merge Process 


Step /.—The initial identification 
list is modified by deleting those 
members for which the freedom has 
been constrained, i.e., let 

v for L = 0 


and 


w for R = 0 


then: 



n = gross matrix size 



L 

W 

i 


0 

2 


u V 

3 

M 

w 

4 


e 

5 


V 

6 


u 

7 

R 

6 

8 


V 

9 


u 

10 


n= 10 


Note: In the program, L, M, and R are specified by 1, 2, and 3, respectively. Also w, 6, v, and u 
are denoted by 1 , 2, 3, and 4, respectively. The merge list is packed. The example from above would 
appear: 


Row 1 

1 

1 

2 

1 

2 

3 

1 

4 

4 

2 

1 

5 

2 

2 

6 

2 

3 

7 

2 

4 

8 

3 

2 

9 

3 

3 

10 

3 

4 
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Step 2.— Given elemental matrix S 

Merge S for left part (L, M) and right part (M, R) 

L M 

w 6 y u w Q v u Column codes 


w 

r"” 0 — 

6 

i 

L 


V 


u 

i 

i 

i 

i 

l 

Hj 

1 

1 

L 

w 

1 

| 

9 

1 

M 

V 

1 b 

i 1—1 

u 

k 1 

1 


u 

t 


Row codes 

a. Member’s , row code L, u — »»-k r = 3 

1 col code L, v -►k c (does not exist) 

Member will be deleted. 

b. Member’s / row code M, v — ►kj. = 6 

1 col code M, 0 —►kg = 8 

Then C(6, 8) = C(6, 8) + S(7, 6) 

where C = array containing merged matrix 

Step 3 . -The left edge spring constant is added to the L, L nodal component, and for the right 
edge the R, R is added to. 
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1.2 OVERLAY STRUCTURE 


Main overlay 0.0 
Program name S0352A 
Subroutines PRINTC 

ERROR 

Primary overlay 1.0 

Primary overlay 2.0 

Primary overlay 4.0 

Program name: DATIN 

Subroutines: PAC 

PRERD 

NEXTC 

UN PAC 

PLTDEF 

RDTBLE 

TBPOINT 

INFORM 

ST R MOV 

PLACE 

COMCHK 

j 

Program name: SELECT 

Program name: 
STRAIN 
Subroutines: 
ASTAR 
SIMEQ 

Secondary 
overlay 2.1 

Secondary 
overlay 2.2 

Secondary 
overlay 2.3 

Program name: 
WAVFND 

Program name: 
GALUP 
Subroutines: 
ASTAR 
REDUC3 
EECM 
SIMEQ 

Program name: 
LOADCL 
Subroutines: 
LOADCN 
SORT 
DC 

MERGEC 

STOREC 

CHMDET 

DBLERT 

AROOT 

PLATEC 

PLTCLC 

CALVET 

RGENER 

SOLVEC 

PAC 

DTC 

CDTM 

VIPDR 

ZARK2 

CXINV 
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1.3 BASIC PROGRAM FLOW 


Figures 1 through 6 are a visual elaboration of the basic control flow of the program. The 
overall flow is indicated on the first figure, with elaboration of each part except GALUP and 
STRAIN on the following flow charts. 
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FIGURE 1. -OVER ALL FLOW CHART 













FIGURE 2. -FLOW CHART OF PRIMARY OVER LA Y 1.0 











FIGURE 3 . -FLOW CHA RT OF SECOND A RY O VERLA Y 2. 1 
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FIGURE 4.-DETAILED FLOW FOR CONFIGURATION 2 IN WA VFND 
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Starting with guesses, loop until roots found 
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FIGURE 6.— DETAILED FLOW CHART OF DC 










1.4 COMMON BLOCK DEFINITIONS 


All common block variables are here defined. Any variables not in common are defined, where 
desirable, by comment cards in the subroutine. 


Common 

block 


Description 

ABD 

A 

Extensional stiffness matrix 


B 

Coupling stiffness matrix 


D 

Bending stiffness matrix 


MATXB 

B matrix switch 

= 0 flat plate has zero B matrix 

= 1 flat plate has nonzero B matrix or curved plate 

ACON 

ACON 

A constant that is set in the DATIN subroutine according to the load option 
and is used in GALUP subroutine for the upper bound calculation 

ANGCM 

NOANG 

Number of angles in the angle table (an input) 


ANGLE 

Angles for fiber direction (angle table) 

BOUND 

SUP 

Upper bound for critical load first as input or as calculated by Galerkin 
method; second, as end result of LOADCN iterations 


SLW 

Lower bound for critical load first as input ; second, as end result of 
LOADCN iterations 


UPBND 

Upper bounds for critical loads as calculated in GALUP by the Galerkin 
method 
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Common 

block 

Variable 

Description h 

CNTROL 

DONSWI 

Donnell switch j 


1PC 

Array of intermediate switches 


JPC 

Array of control switches 

CONPNT 

1PRNT1 

Applied inplane loads check print switch 


1PRNT2 

R matrix check print switch 


IPRNT3 

Xj, Xf and stiffness matrix check print switch 


1PRNT4 

Merged stiffness check print switch 


1PRNT5 

Load search check print switch 

CONTRL 

DOUBLE 

Double root detection switch 


NZERO 

Nonzero determinant switch 


POSTVE 

Positive definite determinant switch 

CWAVAL 

RMA 

Current wave number being searched 

DBTRAN 

LRTRN 

Return address for double root detection 

DCARD 

KTITL 

Title for problem which is placed at top of each output page J 




| Common 
1 block 

— 

Variable 

Description 


1KARD 

Array of last 78 characters on an input card 


1ERFLG 

Input error flag (or switch) -controls flow in program’s input section 


1CRD 

First two characters on an input card 


NPAGE 

Page number for the output of a problem 

DSTIFF 

DSTIFF 

Nodal stiffnesses for plate edges 

ICK 

IX \ 

Count of number of determinant evaluations 

INFO 

WIDE 

Width of flat plate; chord width of a curved plate 


XL1 

First load value (input) 


i 

XL2 

Second load value (input) 


SWAVE 

Half-wave length 


IPRLFT 

Left boundary option (SS, CL, FR, blank) 


IPRRHT 

| Right boundary option (SS, CL, FR, blank) 


KEND 

i 

Pointer to VMIN’s position 


VMIN 

Local minimum buckling load 


1 ID1SP 

Left and right boundary conditions 

LOADVL 

I AN1! 

I Applied inplane load in the 1 1 direction 
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AN 1 2 I Applied inplane shear in the 1 2 plane 


AN22 I Applied inplane load in the 22 direction 


MERGE I NRWS I Number of rows in merged stiffness matrix 


NCLS I Number of columns in merged stiffness matrix 


MRLIST I List of merge information 


NDIR I NDIR I Switch controlling direction of calculations in the WA VEND subroutine 


NND I NND 


Number of negatives on diagonal of merged stiffness matrix 


OPTION I LDOPT I Load search o ption 


BNDLWR f Input lower bound for load search 


BNDUPR ! Input upper bound for load search 


Constant NX load 


Constant NY load 


RIXY Constant NXY load 


NX/NY ratio 


RXXY 1 NX/NXY ratio 


RYXY NY/NXY ratio 
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Variable 


PLATE 

IETYPE 

Element type 
= 1 flat plate 
= 2 curved plate 


PLONG 

Length of plate 


WIDTH 

Plate half width 


RADIUS 

Plate radius 

PRECOM 

KREAD 

Switch signalling whether or not to read a card on initial entry to 
PRERD 


TCARD 

Array of input card images 


I EOF 

Switch indicating whether or not an EOF has been read 

PROOTC 

NR 

Number of a roots to be obtained as a group 


NT 

Total number of a roots obtained 


KXZ 

Switch used in AROOT to indicate in B = 0 case if roots are from 2 x 2 ( = 1 ) 
or 1 x 1 (=0) 


KXT 

Switch used in PLATEC to indicate in B = 0 case if columns arc first four 
(=0) or last four (similar to KXZ (=1 ) 


KXY 

Indicates if B = 0 and flat plate (=1 ) and if not (=0) 


PIX 

Array containing a root type information 

| 

P 

Complex array containing a roots 
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Common 

block 

Variable 

Description ■ 

RC 

R 

Array containing stability equations matrix 

RDCOM | 

CARD § 

: 

Array containing ordered input cards 


KDNT 

Total number of data cards 


KDNO 

Counter for data card image to be processed 

STRCM 

NLAY 

Number of layers 


THTA 

Angle for fiber direction for each layer 


MTOPT 

Type of material properties input 
= 0 Engineering constants 

= 1 Matrices A, B (optional), and D are supplied 

TABLE 

NOTHKV 

Number of items in THICK array 


MTBCNT 

Number of rows in ITABLE array 


THICK 

Lamina thicknesses 


ITABLE 

Plate data control table 

(1.1) Table number 

(1.2) -Number of items in table 


T 

Plate data 

(! ,J) -E-modulus for direction 1 
(2,J)- E-modulus for direction 2 
(3, J) Poisson’s ratio 
(4,J)--G-modulus (G! 2) 

(5,J) fiber angle (degrees) 



Common 

block 

Variable 

Description 

WAVE 

MSOPT 

List search option 

= 1 stop search when local minimum is first detected 
~ 2 search complete list for minimum 


NOWAVE 

Number of wave values in list 


WAVE 

List of wave values to be searched 


IWTYPE 

Type of wave list 
= 1 input as m values 
= 2 input as X values 


MUST 

List of m numbers 


VLOAD 

Array of variable loads obtained for each wave value 

WVSRCH 

NTRY 

Switch for number of critical wave lengths to seek 

0 = only one to find 

1 = seeking first critical with two to find 

2 = seeking second critical with two to find 


NFL1P 

Switch for flipping critical wave length search between two sections 


NHALF 

Section in which wave search being made 


BORDR1 

The lower bound (half-wave length inverse) for critical wave search 


BORDR2 

The upper bound (halt-wave length inverse) for critical wave search 


BORDR3 

The bound beyond which no sharp dropoff in critical load is permitted 


WAVEP 

The critical half-wave length inverse found for first try 


WAVES 

A half-wave length inverse to be used initially on second try 


20 





Co 
block 



A half-wave length inverse to be used initially on second try along with WAVES 


The critical minimum load found for first try 


Critical load value corresponding to WAVES 


Critical load value corresponding to WAVHT 


Critical load corresponding to IIOLDW 


Half-wave length inverse value being held from previous iteration 


ZCOM 


Counter (initially) and length (finally) for arrays of wave lengths, upper bounds, 
and critical loads, which are saved for an output summary 


Half-wave lengths array 


Critical load array 


ZERROR 


Upper bound array 

Return address after error detection during problem execution 



21 







1.5 BASIC DATA FLOW 


Variable 

Common block 

Where set 

Where used 

A 

ABD 

PLTDEF 

GALUP 

PLTCLC 

RGENER 

STRAIN 

ACON 

ACON 

DATIN 

GALUP 

ANGLE 

ANGCM 

DATIN 

PLTDEF 

AN1 1 

LOADVL 

DC 

AROOT 

PLTCLC 

RGENER 

AN 12 

LOADVL 

DC 

AROOT 

PLTCLC 

RGENER 

AN22 

LOADVL 

DC 

DC 

AROOT 

PLTCLC 

RGENER 

B 

ABD 

PLTDEF 

GALUP 

PLTCLC 

RGENER 

STRAIN 

BNDLWR 

OPTION 

DATIN 

LOADCL 

BNDUPR 

OPTION 

DATIN 

SELECT 

LOADCL 

BORDR1 

WVSRCH 

WAVFND 

SELECT 




WAVFND 
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Variable 

Common block 

Where set 

Where used 

B0RDR2 

WVSRCH 

WAVFND 

SELECT 

WAVFND 

B0RDR3 

WVSRCH 

WAVFND 

SELECT 

WAVFND 

CARD 

RDCOM 

PRERD 

PRERD 

DATIN 

NEXTC 

PLTDEF 

RDTBLE 

D 

ABD 

PLTDEF 

GALUP 

PLTCLC 

RGENER 

STRAIN 

DONSWI 

CNTROL 

DATIN 

PLTCLC 




RGENER 

DOUBLE 

CONTRL 

DBLERT 

LOADCN 

DSTIFF 

DSTIFF 

DATIN 

INFORM 

MERGEC 

HOLDL 

WVSRCH 

WAVFND 

WAVFND 

HOLDW 

WVSRCH 

WAVFND 

WAVFND 

ICRD 

DCARD 

DATIN 

DATIN 

IDISP 

INFO 

DATIN 

INFORM 

IEOF 

PRECOM 

S0352A 

PRERD 



STRAIN 

STRAIN 
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Variable 

Common block 

Where set 

Where used 

IERFLG 

DCARD 

DATIN 

DATIN 

IETYPE 

PLATE 

DATIN 

DATIN 

PLTDEF 

INFORM 

IKARD 

DCARD 

DATIN 

PRERD 

PRERD 

1PC 

CNTROL 

DATIN 

LOADCL 

LOADCN 

DC 

IPRLFT 

INFO 

DATIN 

DATIN 

INFORM 

IPRNT1 

CONPNT 

DC 

AROOT 



AROOT 


IPRNT2 

CONPNT 

DC 

RGENER 

IPRNT3 

CONPNT 

DC 

PLATEC 

PLATEC 

IPRNT4 

CONPNT 

LOADCN 

DC 

IPRNT5 

CONPNT 

LOADCL 

LOADCN 

IPRRHT 

INFO 

DATIN 

DATIN 

INFORM 
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Variable 

Common block 

Where set 

Where used 

ITABLE 

TABLE 

DATIN 

PLTDEF 

IWTYPE 

WAVE 

DATIN 

DATIN 

INFORM 

LOADCL 

IX 

ICK 

LOADCL 

LOADCN 



DC 

DC 

JPC 

CNTROL 

DATIN 

S0352A 

DATIN 

SELECT 

KDNO 

RDCOM 

DATIN 

DATIN 



NEXTC 

NEXTC 



PLTDEF 

PLTDEF 



RDTBLE 

RDTBLE 

KDNT 

RDCOM 

PRERD 

PLTDEF 

RDTBLE 

KEND 

INFO 

WAVFND 

LOADCL 

LOADCL 

KREAD 

PRECOM 

S0352A 

PRERD 

PRERD 

KTITL 

DCARD 

PRERD 

PLTDEF 




LOADCL 

KXT 

PROOTC 

PLATEC 

PLTCLC 

KXY 

PROOTC 

AROOT 

AROOT 



PLTCLC 

PLTCLC 

RGENER 

DTC 
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Variable 

Common block 

Where set 

Where used 

KXZ 

PRODTC 

AROOT 

AROOT 

RGENER 

DTC 

LDOPT 

OPTION 

DATIN 

DATIN 

INFORM 

GALUP 

DC 

STRAIN 

LRTRN 

DBTRAN 

DG 

DBLERT 

LTRAN 

ZERROR 

S0352A 

ERROR 

MATXB 

ABD 

PLTDEF 

AROOT 

MLIST 

WAVE 

DATIN 

DATIN 

INFORM 

LOADCL 

MRLIST 

MERGE 

DATIN 

MERGEC 

MSOPT 

WAVE 

DATIN 

■/ 

DATIN 

SELECT 

LOADCL 

MTBCNT 

TABLE 

DATIN 

PLTDEF 

MTOPT 

STRCM 

PLTDEF 

STRAIN 

NCL.S 

MERGE 

DATIN 

LOADCL 

MERGEC 

NDIR 

NDIR 

SELECT 

WAVFND 

SELECT 

WAVFND 

LOADCL 
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Variable 

Common block 

Where set 

Where used 

NFLIP 

WVSRCH 

WAVFND 

WAVFND 

NHALF 

WVSRCH 

WAVFND 

WAVFND 

NLAY 

STRCM 

PLTDEF 

STRAIN 

NLCNT 

ZCOM 

SELECT 

LOADCL 

LOADCL 

NND 

NND 

DC 

DC 

NOANG 

ANGCM 

DATIN 

PLTDEF 

NOTHKV 

TABLE 

DATIN 

PLTDEF 

NOWAVE 

WAVE 

DATIN 

DATIN 



WAVFND 

SELECT 

GALUP 

LOADCL 

NPAGE 

DCARD 

PRERD 

PRERD 



PLTDEF 

PLTDEF 



LOADCL 

LOADCL 

NR 

PROOTC 

DC 

AROOT 

AROOT 

NRWS 

MERGE 

DATIN 

LOADCL 

MERGEC 

NT 

PROOTC 

AROOT 

AROOT 

NTRY 

WVSRCH 

WAVFND 

WAVFND 

NZERO 

CONTRL 

DC 

LOADCN 

P 

PROOTC 

AROOT 

AROOT 

PLATEC 


27 



Variable 

Common Hock 

Where set 

Where used 

FIX 

PRQGIC 

AROOT 

AROOT 

PLONG 

PLATE 

DATIN 

INFORM 

POSTVE 

CONTRL 

DC 

LOADCN 

R 

RC 

RGENER 

RGENER 

DTC 

RADIUS 

PLATE 

DATIN 

DATIN 

INFORM 

GALUP 

PLTCLC 

RGENER 

RIX 

OPTION 

DATIN 

INFORM 

DC 

RIXY 

OPTION 

DATIN 

INFORM 

DC 

RIY 

OPTION 

DATIN 

INFORM 

DC 

RMA 

CWAVAL 

LOADCL 

PLTCLC 

RGENER 

RXXY 

OPTION 

DATIN 

INFORM 

DC 

STRAIN 

RXY 

OPTION 

DATIN 

INFORM 

DC 

STRAIN 

RYXY 

OPTION 

DATIN 

INFORM 

DC 

STRAIN 

SLW 

BOUND 

LOADCL 

LOADCL 
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Variable 

Common Mock 

Where set 

Where used 

SUP 

BOUND 

LOADCL 

LOADCL 

SWAVE 

INFO 

DATIN 

INFORM 

LOADCL 

T 

TABLE 

DATIN 

PLTDEF 

TCARD 

PRECOM 

PRERD 

PRERD 

THICK 

TABLE 

DATIN 

PLTDEF 

THTA 

STRCM 

PLTDEF 

STRAIN 

TLOAD 

ZOOM 

LOADCL 

LOADCL 

TWAVE 

ZOOM 

LOADCL 

LOADCL 

UPBND 

BOUND 

GALUP 

LOADCL 

UPPER 

ZOOM 

LOADCL 

LOADCL 

VLOAD 

WAVE 

WAVFND 

WAVFND 



LOADCL 

LOADCL 

VLOADP 

WVSRCH 

WAVFND 

WAVFND 

VLOADS 

WVSRCH 

WAVFND 

WAVFND 

VLOADT 

WVSRCH 

WAVFND 

WAVFND 

VMIN 

INFO 

LOADCL 

LOADCL 

STRAIN 

WAVE 

WAVE 

DATIN 

DATIN 



SELECT 

SELECT 



WAVFND 

WAVFND 

GALUP 
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Variable 

Common block 

Where set 

Where used 

WAVEP 

WVSRCH 

WAVFND 

WAVFND 

WAVES 

WVSRCH 

WAVFND 

WAVFND 

WAVET 

WVSRCH 

WAVFND 

WAVFND 

WIDE 

INFO 

DATIN 

INFORM 

WIDTH 

PLATE 

DATIN 

DATIN 

INFORM 

GALUP 

XL1 

INFO 

DATIN 

DATIN 

XL2 

INFO 

DATIN 

DATIN 
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1.6 PROGRAM AND SUBPROGRAM DESCRIPTIONS 


Each subroutine is described in this section. The descriptions are placed in alphabetical 
order. 

Contributions made by Paul Lu, Claude Gagnon, and Edith Shook of Boeing Computer 
Services are gratefully acknowledged. 

Subroutine AROOT 

Purpose: To find the roots of the characteristic equations associated with the stability 

equation (Eq. A.37, ref. 1) 

Called From: DC 

Input: (a) Calling sequence— none 

(b) COMMON-/ABD/, MATX B 

/CONPNT/, IPRNT1 
/LOADVL/ AN1 1, AN12, AN22 
/PROOTC/ NR 

Discussion: The characteristic roots are obtained by a root finder. The roots are then 

examined for conjugate pairs and the occurrence of double roots. For flat 
plates with zero B matrix the problem is uncoupled into two fourth-order 
problems. 

Routined Called: ZARK2 

DBLERT 
ERROR 

Output: (a) Calling sequence— none 

(b) COMMON— /PROOTC/, NR, NT, KXZ, KXY, PIX, P 

(c) User I/O files output-check print of load values and characteristic roots. 

Error message on root searching. 
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Subroutine ASTAR (ASTAR (AST, A, B, D, IBAD) 


Purpose: 

Called From: 

Input: 

Discussion: 
Routines Called: 

Output: 


To perform calculations AST = A - B*D‘* * B 

STRAIN 

GALUP 

Calling sequence 

A— extensional stiffness matrix (3 x 3) 
B— coupling stiffness matrix (3 x 3) 

D— bending stiffness matrix (3 x 3) 

None 

SIMEQ 


AST— reduced matrix 

IBAD— error return (0 = ok, 1 = error) 
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Subroutine CAL VET (A, ER, IPR) 


Purpose: 

Called From: 
Input: 

Discussion: 

Routines Called: 
Output: 


To solve a 3 x 3 homogeneous complex system of equations A * X = 0 by 
inverse iteration technique. 

PLTCLC 

Calling sequence 

A— 3 x 3 complex array containing given matrix 

The given square matrix A is decomposed into lower and upper triangular 
matrices, L and U, by Crout’s method with partial pivoting. Then Wielandt’s 
inverse iteration for finding the eigenvector is used to generate the required 
solution vector which is normalized to the third component. 

None 

Calling sequence 

V- complex array returning the solution vector 
IPR- working array for pivoting information 
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Function CDTM (CA, NR, N, V) 


Purpose: 
Called From: 
Input: 


Discussion: 


Routines Called: 
Output: 


To evaluate the determinant of a complex square matrix 
DTC 

(a) Calling sequence 

A— elements of matrix stored in a complex array 

NR— 2 x NC (NC-maximum row dimension of complex array A) 

N— dimension of the square matrix 
V— a scratch array of length > N 

(b) COMMON-none 

The given square matrix A is decomposed into lower and upper triangular 
matrices, L and U, by Crout’s method with partial pivoting and row equili- 
bration. So we have n 

det (PA) - det (L) det (U) = ?r ®i,i 

i = 1 

where PA = LU and P is a product of permutation matrices, and 

det (A) = (-l) k S 
i = 1 

where k is the total number of row permutations performed on A. The 
routine uses the standard FORTRAN convention for storing complex matrices, 
but it does not use FORTRAN complex arithmetic. The magnitude of real 
and imaginary part of the determinant value must be between the lower and 
the upper bounds of the floating point numbers on the machine. CDTM is 
typed complex. 

VIPDR 

Y = CDTM (CA, NR, N, V) 

Y— the complex determinant value 
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Subroutine CHMDET (A, NR, N, DN, D, ID, NND) 


Purpose: 

Called From: 
Input: 

Discussion: 


Routines Called: 


To compute the determinant of a complex Hermitian matrix, and to find the 
number of negative diagonal elements. 

DC 

Calling sequence 

A— two-dimensional complex array (The upper triangular portion contains 
the upper triangle of the given Hermitian matrix.) 

NR— maximum row dimension of the complex array A 
N— actual order of the given Hermitian matrix 

A complex analog of the symmetric gaussian elimination algorithm is used. 

The given Hermitian matrix A is decomposed to produce the diagonal ele- 
ments djj of the matrix D as follows: 


A = LDU (1) 

where D is a diagonal matrix, L is a unit complex lower triangular matrix, 
and U = L* is the conjugate transpose of L. The diagonal elements d^ are 
real because the complex matrix A is Hermitian. 

Since det (L) n = det (U) - 1 , 

det (A) = 7 r d^ (2) 

i = 1 

In the routine CHMDET the product in equation (2) is expressed in a form 
det (A) = a • 2^ where 1/16 < | a | < 1 , to avoid underflow or overflow. 
The number of negative diagonal elements is counted in the process of 
decomposition. 

None 
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Output: 


A— the lower triangular portion contains the matrix L without diagonal and 
the upper triangular portion is left unchanged 
DC— the reciprocals of the elements of the diagonal matrix D 

D— the fractional part a of the determinant - a • 2 . It is set equal to zero 

if decomposition failed. 

ID— the exponential part b of the determinant - a • 2^ 

NND— number of negative diagonal elements 
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Subroutine COMCHK (CARD, NCOM) 


Purpose: 
Called from: 
Input: 

Discussion: 

Output: 


To check certain input cards for commas and slashes 
PRERD 

Calling sequence 

CARD— array of Hollerith data from a data input card 

A data input card that is normally fixed field is checked for commas 
and slashes to determine if it has been input as a free field card. 

Calling sequence 

NCOM— an indicator whose value is 0 (no commas or slashes) or 1 
(commas or slashes). 
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Program DATIN 


Purpose: 
Called From: 
Input: 

Discussion: 


Routines Called: 
Output: 


To read a data set and set up the problem from it to be worked 
S0352A 

(a) Calling sequence— none 

(b) COMMON — /RDCOM/ CARD, KDNT 

The following functions are performed on the preread cards: 

(1) C2 card content is placed into IPC array 

(2) C3 card content is placed into JPC array 

(3) Load case information constructed LC card 

(4) W 1 , W2, W3 cards supply wave search information 

(5) BL and BR supply the information for construction of merge list 

(6) Tables for lamina thickness, material properties, and fiber angles are 
processed 

(7) PI card is processed 

PRERD 

SECOND 

RDTBLE 

ASIN 

PLTDEF 

(a) Calling sequence— none 

(b) COMMON— /ANGCM/, NOANG, ANGLE 

/CNTROL/, DONSWI, IPC, JPC 
/DCARD/, IKARD, IERFLG 

/INFO/, WIDE, XL1 , XL2, SWAVE, IPRLFT, IPRRHT 
/PLATE/, lETYPE, PLONG, WIDTH, RADIUS 
/TABLE/, NOTHKV, MTRCNT, THICK, IT ABLE, T 
/DSTIFF/, DSTIFF 
/MERGE/, NRWS, NCLS, MRLIST 

/OPTION/, LDOPT, BNDLWR, BNDUPR, RIX, RIY, RIXY, 
RXY, RXXY, RYXY 


38 



/WAVE/, MSQPT, NOWAVE, WAVE, IWTYPE, MUST 
/AGON/, AGON 
/RDCOM/, KDNO 

(c) User I/O files output— diagnostics and error messages 
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Subroutine DBLERT 


Purpose: 

Called From: 
Input: 

Discussion: 

Routines Called: 
Output: 


To serve as error trap for double root detection 
AROOT 

(a) Calling sequence— none 

(b) COMMON— /DBTRAM/LRTRN 

In a double root situation, the condition is flagged and return is made to load 
control routine to attempt to avoid the problem. 

None 

(a) Calling sequence— none 

(b) COMMON — /CONTRL/ , DOUBLE 

(c) User I/O files output— error message indicated double root condition. 
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Function DC (VLOAD, DBMA, Ml, Nl, DX, IDX) 


Purpose: 

Called From: 
Input: 


Discussion: 
Routines Called: 

Output: 


To calculate the load configuration for a given load option and varying load 
component and to control the formation of and the determinant calculation 
of the resulting stiffness matrix. 

LOADCL through LOADCN 

(a) Calling sequence 

VLOAD— varying load component 
DBMA— work array for determinant calculation 
Ml— row dimension of DBMA array 
Nl —column dimension of DBMA array 

(b) COMMON— /CNTROL/, IPC 

/CONPNT/, IPRNT4 
/CONTRL/, DOUBLE 
/ICK/, IX 

/OPTION/, RIX, RIY, RIXY, RXY, RXXY, RYXY 


The subprogram determines conditions that might occur in the determinant 
calculation such as: zero determinant and double roots among characteristic 
equation associated with stability equations. 

PLATEC 

MERGEC 

PRINTC 

CHMDET 

(a) Calling sequence 

DX— the exponential part b of the determinant a • 2^ 

IDX— the fractional part a of the determinant a * 2^ 

(b) COMMON- /CONPNT/, IPRNT1, IPRNT2, IPRNT3 

/CONTRL/, NZERO, POSTVE 
/DBTRAN '/, LRTRN 
/ICK/, IX 

/LOADVL/, AN 11, AN12, AN22 
/PROOTC/, NR 
/NND/, NND 

(c) User I/O files output— summary of determinant calculation 
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Function DTC (P) 


Purpose: 

To obtain, for a complex number P, the determinant of the matrix associated 
with the stability equations (eq A.37, ref. 1). 

Called From: 

AROOT 

Input: 

(a) Calling sequence— P, a complex number 

(b) COMMON — /PROOTC/ , KXY, KXZ 

IRC/, R 

Discussion: 

For a complex number P, the matrix for the stability equations is generated 
by RGENER. Next, depending on the element type, the determinant for a 
3 x 3 or 2 x 2 or 1 x 1 matrix is obtained. 

Routines Called: 

CDTM 

RGENER 

Output: 

Y = DTC (P) 

Y— value of determinant 
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Subroutine ERROR (NAME, NO) 


Purpose: 

Called From: 
Input: 

Discussion: 
Routines Called: 
Output: 


To handle fatal error conditions 

SOLVEC, LOADCN, STRAIN, AROOT, GALUP 

Calling sequence 

NAME— name of routine where error was detected 
NO— error number 

The routine name and error number is printed. A return to S0352A is made. 
None 

(a) Calling sequence— none 

(b) COMMON— none 

(c) User I/O files ouput— routine name and error number 


43 



Program GALUP 


Purpose: 

Called From: 
Input: 


Discussion: 


Routines Called: 
Output: 


To calculate an upper bound for a clamped plate by the Galerkin method. This 
upper bound is used in the calculation procedure which finds the critical load. 

SELECT 

(a) Calling sequence— none 

(b) COMMON- /ABD/, A, B, D 

/PLATE/, WIDTH, RADIUS 
/WAVE/, NOWAVE, WAVE 
/ACON/, ACON 
/OPTION/, LDOPT 


The A, B, and D, matrices are used for the calculation of required constants. If 
the Bj g term or B2 g term is nonzero, a modified B matrix is determined for 
use by GALUP. Required constants differ according to the load option in use. 

In essence the Galerkin method results in the following matrix equation: 

(A + X B) j W j =0 

The minimum positive eigenvalue is found. This is used for the upper bound. 

ASTAR 

REDUC3 

EECM 

GALUP 

(a) Calling sequence— none 

(b) COMMON— /BOUND/, UPEND 
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Subroutine INFORM 


Purpose: To summarize the input options. 

Called From: PLTDEF 

Input: (a) Calling sequence— none 

(b) COMMON— /DCARD/, KTITL, NPAGE 
/DSTIFF/, DSTIFF 

/OPTION/, LDOPT, BNDLWR, BNDUPR, RIX, RIY, RIXY 
RXY, RXXY, RYXY 

/PLATE/, IETYPE, PLONG, WIDTH, RADIUS 

/WAVE/, IWTYPE, MLIST 

/INFO/, WIDE, SWAVE, IPRLFT, IPRRHT 


Discussion: 

None 

Routines Called: 

None 

Output: 

(a) Calling sequence— none 


(b) COMMON-none 


(c) User I/O files output— summary print 


45 



Program LOADCL 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


To control the search for a minimum buckling load value over a range of 
half-wave values. 


' SELECT 

(a) Calling sequence— none 

(b) COMMON-/BOUND/, UPBND 

/CNTROL/, IPC 

/MERGE/, NRWS, NCLS 

/OPTION/, BNDLWR, BNDUPR 

/WAVE/, MSOPT, NOWAVE, WAVE, IWTYPE, MLIST 

/DCARD/, KTITL, NPAGE 

/NDIR/, NDIR 


The search runs on basis of saying either the minimum is first half-wave value 
for which the next wave value has a larger buckling load, or the minimum is 
the half-wave value with the smallest buckling load where all wave values are 
searched. 

LOADCN 

DC 

(a) Calling sequence— none 

(b) COMMON —/BOUND/, SUP, SLW 

/CONPNT/, IPRNT5 

/CWAVAL/, RMA 

/ICK/, IX 

/WAVE/, VLOAD 

/INFO/, SWAVE, KEND, VMIN 

(c) User I/O files output— summary of minimum load search 
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Subroutine LOADCN (DB, DBMA, Ml, Nl, IP, SUP, SLW) 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


For a given interval, to find the buckling load for the interval to a specified 
tolerance. 

LOADCL 

(a) Calling sequence 

DB— name of function program which calculates determinant 

DBMA— work array for determinant 

Ml— row size of determinant 

N 1 —column size of determinant 

SUP— upper end of given interval 

SLW— lower end of given interval 

(b) COMMON— /CNTROL/, IPC 

/CONPNT/,IPRNT5 

/CONTRL/, DOUBLE, NZERO, POSTVE 
/ICK/, IX 

The load search uses a modified bisection technique. When a 0-1 root-crossing 
situation is detected, an effort is made to improve convergence by choosing 
a trial load which is nearer the actual buckling load. The program also handles 
abnormal conditions such as double roots and zero determinant by perturbing 
the trial load. 

ERROR 

(a) Calling sequence 

SUP— upper value of calculated interval containing buckling load 
SLW— lower value of calculated interval containing buckling load 

(b) COMMON — /CONPNT / , IPRNT4 

(c) User I/O files output— check print of load search 
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Subroutine MERGEC (A, M, S, IP) 


Purpose: 
Called From: 
Input: 


Discussion: 


Routines Called: 
Output: 


To merge an elemental stiffness matrix into the global stiffness matrix. 

DC 

(a) Calling sequence 

A— complex array containing global stiffness matrix 
M— row dimension of array A 

S— complex array containing elemental stiffness matrix 
IP— pointer to which side elemental matrix is on. 

= 1 Left side 
2 Right side 

(b) COMMON— /DSTIFF/, D STIFF 

The elemental matrix is merged in by nodal components with the diagonal 
component being modified by the appropriate nodal stiffness matrix (if any). 




— — Nodal component 






Elemental matrix (8 x 8) 


STOREC 
Calling sequence 

A— complex array containing global stiffness matrix 
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Function NEXTC (IN, IGUT, IA, IPOS) 


Purpose: 

Called From: 
Input: 

Discussion: 
Routines Called: 
Output 


To read table pointer card as needed and return one character each time entered. 
TBPOINT 

(a) Calling sequence 

IPOS— location of character in string 

(b) COMMON— /RDCOM/, CARD, KDNO 
IPOS must be >78 for first entry. 

None 

(a) Calling sequence 

NEXTC— character status flag 

if IA = number, NEXTC = 1 R9 
- character, NEXTC = IA 
IA— the IPOS character from the string 

(b) COMMON-/RDCOM/, KDNO 


49 



Subroutine PAC (A, I, J, B) 


Purpose: 

Called From: 

Input: 

Discussion: 
Routines Called: 
Output: 


To pack the (J-I+l ) right-most bits of the word B into bit positions 1-3 of the 
word A with the remaining bits unchanged. 

DATIN 

STRMOV 

STOREC 

Calling sequence 
I— left position 
J— right position 

B— word from which packing is being done 
Routine is in COMPASS using RUN linkage conventions. 

None 

Calling sequence 

A— word into which packing is being done 
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Subroutine PLACE (A, DUM) 


Purpose: 
Called from: 
Input: 

Discussion: 

Output: 


To place input data into A, B, or D array. 


PLTDEF 


Calling sequence 

DUM— array of input values 


None 


Calling sequence 

A— array filled with input values 
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Subroutine PLATEC (ST, XD, XF) 


Purpose: 

Called From: 
Input: 

Discussion: 

Routines Called: 
Output: 


To control the generation of columns of the buckling displacements and 
corresponding loads and produce the resulting stiffness matrix. 

DC 

(a) Calling sequence— none 

(b) COMMON-/CONPNT/, IPRNT3 

/PROOTC /, NT, KXY 

The equations for buckling displacements as defined by equation A.45 (ref. 1) 
and corresponding loads defined by equation A.5 1 (ref. 1 ) are generated and 
the stiffness matrix is obtained by the relationship specified in equations A.5 2 
and A.5 3 (ref. 1). 

PLTCLC 

PRINTC 

SOLVEC 

(a) Calling sequence 

ST— 8 x 8 complex array containing stiffness matrix 
XD— 8 x 8 complex array containing buckling displacements 
XF-8 x 8 complex array containing corresponding loads 

(b) COMMON-/PROOTC/, KXT 
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Subroutine PLTCLC (P, DL, FL, DR, FR) 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


For a P value, to generate the columns of the deflection and force matrix for 
the right and left sides of the plate. 

PLATEC 

(a) Calling sequence— P - P value 

(b) COMMON-/ ABD/, A, B, D 

/CNTROL/, DONSWI 
/CWAVAL/, RMA 
/LOADVL/, AN1 1, AN12, AN22 
/PLATE/, WIDTH, RADIUS 
/RC/, R 

/PROOTC/, KXT, KXY 

The deflection matrix represents equations A.40, A.41 , A.42 and A.43 of 
reference 1. The force matrix represents equations A.46, A.47, A.48, and A.49 
of reference 1 . In the case of B = 0 and flat plate the equations have been modi- 
fied accordingly. 

RGENER 

CALVET 

(a) Calling sequence 

DL— complex matrix containing left-side deflection column 
FL— complex matrix containing left-side force column 
DR— complex matrix containing right-side deflection column 
FR— complex matrix containing right-side force column 

(b) COMMON —/PROOTC/ , KXY 
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Subroutine PLTDEF (MATOPT, INOPT, LAYERS, SM2X, ZNX, IER) 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


To produce the plate element extensional, coupling, and bending stiffness 
matrices. 

DATIN 

(a) Calling sequence 

MATOPT— type of material properties input 
INOPT— input mode when MATOPT = 0 
LAYERS— number of laminas in plate 

(b) COMMON— /PLATE/, IETYPE 

/TABLE/, NOTHKV, MTBCNT, THICK, ITABLE, T 

/DCARD/, KTITL, NPAGE 

/RDCOM/, CARD, KDNT, KDNO 

/ANGCM/, NOANG, ANGLE 

/CNTROL/, IPC 

The A, B, D matrices are either input directly or calculated from material 
properties input by P2, P3, or LM cards using equations A.4, A.8, A.9, and 
A. 10 (ref. 1). 

TBPOINT 

INFORM 

PLACE 

(a) Calling sequence 

IER— status flag (0 = no error) ( 1 = errors) 

(b) COMMON— /ABD/, A, B, D, MATXB 

/STRCM/, NLAY, THTA, MTOPT 

(c) User I/O files output— summary of element input, print of A, B, and D 

matrices 
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Program PRERD 


Purpose: 

Called From : 
Input: 

Discussion: 
Routines Called: 
Output: 


To provide a preliminary processing of the input cards. 

DATIN 

(a) Calling sequence— none 

(b) COMMON— /RDCOM/ KREAD, TCARD, IEOF 

PRERD brings in a data set from input cards, reorders the data set to a 
specified order, and checks to see if all required input cards were present. 

DATE 

COMCHK 

(a) Calling sequence— none 

(b) COMMON— /DCARD/ KTITL 

/RDCOM/ KREAD, IEOF 
/PRECOM/ CARD, KDNT, KDNO 


55 



Subroutine PRINTC (A, NROW, NCOL, TITLE, NIZ) 


Purpose: 

Called From: 

Input: 

Discussion: 
Routines called: 
Output: 


To print a rectangular complex matrix stored in the array A. 
DC 

PLATEC 
Calling sequence 

A— two-dimensional complex array containing matrix 
NROW— number of rows in matrix 
NCOL— number of columns in matrix 
TITLE— one word containing matrix title 
NIZ— maximum row dimension 


The matrix is printed row wise, four pairs per row. 
None 


(a) Calling sequence— none 

(b) COMMON-none 

(c) User I/O files output— print of matrix 
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Subroutine RDTBLE (IU, ITBL, CARD, T, MAX, L, IER) 


Purpose: 

Called From: 
Input: 

Discussion: 

Routines Called: 
Output: 


To read floating point or integer data from columns 6-75 of card format . . . , 
array with format value, . . ./ 

DATIN 

(a) Calling sequence 

KARD -array containing first card image 
MAX-maximum size of array T 
ITBL— name of Table T 

(b) COMMON-/RDCOM/, CARD, KDNQ, KDNT 

Delimiters are commas, end of physical card, and slash (/). The number of 
significant digits is limited to 14. Blanks are ignored. The incoming values are 
reformatted so they can be read by a FORTRAN format statement. 

STRMOV 

T— array containing floating point numbers 
L— number of values stored in T 
IER— status flag 
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Subroutine REDUC3 (A, B, D, N, NR, IER) 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


To reduce the generalized complex eigenproblem Ax = ABx, where A is a 
Hermitian matrix and B is a positive definite Hermitian matrix, to a complex 
Hermitian eigenproblem CY - AY 

GALUP 

Calling sequence 

A— a two-dimensional complex array (The upper triangular portion contains 
the upper triangle of the given Hermitian matrix A. If the given matrix 
is real, the imaginary part of A should be set to zero.) 

B— a two-dimensional complex array containing the given matrix B, similar 
to setup for matrix A 

NR— maximum row dimension of the complex arrays A and B 

N— actual order of the given matrix 


First, the complex analog of Cholesky decomposition is applied to B, B = LL* 
where L is a lower triangular matrix and L* is the conjugate transpose of L. By 
forward substitutions, the transpose of the upper triangle of the matrix G is 
formed, LG = A. Finally, the upper triangle of the required complex Hermitian 
matrix is computed by solving the equation CL* = G. 

None 

Calling sequence 

A— the desired complex Hermitian matrix C is in the upper triangular portion 
B— the upper triangular portion is as input and the lower triangular portion 
contains the matrix L without diagonals 
D— the real diagonal elements of the lower triangular matrix L 
IER— status flag 

= 1 successful return 

= 0 the B matrix appears not to be positive definite 
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Subroutine RGENER (ALPHA) 


Purpose: 

To generate coefficients of R matrix (eq. A.37, ref. 1 ) for a complex number 
ALPHA and a load configuration. 

Called From: 

DTC 

PLTCLC 

Input: 

(a) Calling sequence— ALPHA— complex number 

(b) COMMON-/ ABD /, A, B, D 

/CONPNT/, IPRNT2 
/CONTRL/, DONSWI 
/CWAVAL/, RMA 
/LOADVL /, AN 1 1 , AN 1 2, AN22 
/PLATE/, RADIUS 
/PROOTC/, KXZ, KXY 

Discussion: 

Depending on element type, a 3 x 3 matrix is generated (curved plate or plate 
with B ^ 0). If B = 0 and flat plate, a 2 x 2 submatrix is generated in upper 
2 x 2 of R and 1 x 1 submatrix is generated in 3, 3 spot of R. 

Routines Called: 

None 

Output: 

(a) Calling sequence— none 

(b) COMMON— /RC/, R 

(c) User I/O files output— check print of R matrix 
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Program SELECT 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


SELECT is the monitor program which choses the various sections of load 
calculation overlay to be executed. 

S0352A 

(a) Calling sequence— none 

(b) COMMON — /CNTROL/, JPC 

/WAVE/, MSOPT, NOWAVE, WAVE 
/OPTION/, BNDUPR 
/NDIR/, NDIR 

/WVSRCH/, BORDR1, BORDR2 
Detailed functions of load overlay 

(1) Picking wave values for use in finding a minimum critical load 

(2) Upper-bound calculation 

(3) Critical load calculation 

WAVFND 

GALUP 

LOADCL 

(a) Calling sequence— none 

(b) COMMON— /NDIR/, NDIR 

/WAVE/, WAVE 
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Subroutine SOLVEC (X, A, B, N, M) 


Purpose: 

Called From: 
Input: 

Discussion: 
Routines Called: 
Output: 


To solve matrix equation X * A = B. 
PLATEC 


Calling sequence 

A— complex array containing A matrix ( M x M order) 

B— complex array containing B matrix-fN x M order) 

N— number of rows in B matrix, N < 8 
M— number of rows in A matrix, M < 8 

To utilize a standard complex linear equation solver, the original problem is 
transformed to an equivalent problem (A^ X ^ = B^). X^ is solved for and then 
transformed to obtain the result, X. 

PRINTC 

ERROR 

CXINV 

X— complex array containing X matrix (result) 

If an apparent singularity is encountered in the decomposition of A^, a call to 
ERROR is made. 
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Subroutine SORT(TWAVE, UPPER, TLOAD, NLCNT, ALOAD, AUPPER) 


Purpose; 

Called From: 
Input: 

Discussion: 
Routines Called: 
Output: 


To reorder the half-wave lengths and corresponding minimum 
load values in ascending order of the half-wave lengths. 

LOADCL 

Calling sequence : 

TWAVE— array of half-wave lengths 
UPPER— array of upper bounds 
TLOAD— array of critical loads 
NLCNT— size of the arrays 

A simple sort is accomplished 

None 


Calling sequence : 

TWAVE, UPPER, and TLOAD reordered 
ALOAD— minimum critical load 

AUPPER— upper bound associated with minimum critical 
load 
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Subroutine STOREC (A, M, W, I, J) 


Purpose: 

Called From: 
Input: 


Discussion: 

Routines Called: 
Output: 


To merge nodal components of element stiffness matrix into global stiffness 
matrix. 

MERGEC 

(a) Calling sequence 

A— array containing global stiffness matrix 
M— row dimension of global stiffness matrix 
W— array containing nodal component 
I— left node number of nodal component 
J— right node number of nodal component 

(b) COMMON— /MERGE/, NEWS, NCLS, MRLIST 

The nodal component node numbers are looked up in global node list to obtain 
the component merge position. 


PAC 

Calling sequence 

A— array containing global stiffness matrix 


63 



Program STRAIN 


Purpose: 
Called From: 
Input: 


Discussion: 
Routines Called: 

Output: 


To process the calculation of the critical strains. 
S0352A 


(a) Calling sequence— none 

(b) COMMON-/ABD/, A, B, D 

/OPTION/, LDOPT, RIX, RIY, RIXY, RXY, RXXY, RYXY 
/STRCM /, NLAY, THTA, MTOPT 
/INFO/, VMIN 
/PRECOM/, IEOF 


The strains are calculated using equation A.56 (ref. 1) 


ASTAR 

SIMEQ 

ERROR 

(a) Calling sequence— none 

(b) COMMON— /PRECOM/, IEOF 

(c) User I/O files output— strain summary print 
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Subroutine STRMOV (SI, Cl, N, S2, C2) 


Purpose: 
Called From: 
Input: 


Discussion: 
Routines Called: 

Output: 


To move a string of characters from one array to another, 

RDTBLE 
Calling sequence 

S 1 -array containing the string to be moved 

Cl —an integer specifying the position of the first character to be moved 
relative to the first position of the array 
N— number of characters to be moved 

C2— an integer specifying the position of the first character is S2 which is 
to receive the characters from SI . 


None 

PAC 

UNPAC 


Calling sequence 

S2— Receiving array name 
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Program S0352A 


Purpose: 

Called From: 
Input: 
Discussion: 
Routines Called : 

Output: 


S0352A is the monitor propam which chooses the various sections of the 
analysis to be used. 


COMMON/CONTROL/, JPC 
COMMON/PRECOMA IEOF 

The propam terminates by reading an end of file on input. 

DATIN 

SELECT 

STRAIN 

EXIT 

(a) Calling sequence— none 

(b) COMMON-none 
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Subroutine TBPOINT (IN, IOUT, MAX, IPN, IER) 


Purpose: 

Input: 

Discussion: 


Routines Called: 
Output: 


To construct a table of pointers from the P3 or LM card. 

Calling sequence 

MAX— required length of table 

Table of pointers may have the following format: 

(I ] , 1 2 , ’ ’ *3/^ 1 » » 5 2 , ^4^1’ > > > ^2^ 

Between each slash there must be the same number of entries. The result 


IPN (1,1) = I 1 , 

(1,2) = Jj, 

(1,3) = K| 

(2,1) = I 2 , 

(2,2) - Jj, 

(2,3) - Kj 

(3,1) = I 2 , 

(3,2) = J 2 , 

(3,3) = 

(4,1) = I 2 , 

(4,2) = J 3 , 

(4,3) -Kj 

(5,1) = I 3 . 

(5,2) = J 4 , 

(5,3) = K 2 


NEXTC 


(a) Calling sequence 

IPN— array containing table of pointers 

(b) COMMON— name 

(c) User I/O files, output— error messages 



Subroutine UNPAC (A, I, J, B) 


Purpose: 

Called From: 
Input: 

Discussion: 
Routines Called; 
Output: 


To extract bits I-J from word A and place them right adjusted into the word B 
with the rest of B being set to zero. 

STRMOV 

Calling sequence 

A— source of bits to be placed in B 
I— first bit position of string from A 
J— last bit position of string from A 

UNPAC is coded in COMPASS using RUN linkage convention. 

None 


Called sequence 

B— word receiving extracted bit string from A 
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Subroutine VIPDR (A, INCA, B, INCB, N, C) 


Purpose: 
Called From: 
Input: 


Discussion: 

Routines Called : 
Output: 


To calculate the inner product of two complex vectors in double precision. 

CDTM 

(a) Calling sequence 

A— first location of vector 1 

INCA— storage increment between successive elements of vector 1 
B— first location of vector 2 

INCB— storage increment between successive elements of vector 2 
N— number of elements in vector 1 or vector 2 

(b) COMMON— none 

For efficiency the routine is coded in COMPASS using RUN linkage conventions. 
The result is stored in C, C + 1 (C should be declared double precision in the 
calling program). 

None 

(a) Calling sequence 

C— first location of double precision result 

(b) COMMON— none 
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Program WAVFND 


Purpose: The WAVFND subroutine selects the wave values (half-wave length inverses) 

to be used in the search for a local minimum critical load. 


Called From: SELECT 


Input: (a) Calling sequence— none 

(b) COMMON -»/NDIR/ , NDIR 

/WAVE/, WAVE, VLOAD 


Discussion: Given three wave values in the data input, the WAVFND subroutine first 

establishes upper and lower limits on the wave values that it will select. A 
third limit is established. Beyond this third limit a sudden dropoff in any 
critical load value found will terminate the search for a minimum critical load. 

Critical load values are calculated elsewhere in the BUCLAP2 program. On 
return to WAVFND for a second pass, the set of three wave value-critical load 
value pairs is classified as to which of four configurations it belongs. See the 
following diagram. 


Configurations 


V 

| Load 

a 

| Load 

V 

| Load 

A 

| Load 

Wave value 

Wave value 

Wave value 

Wave value 

1 

2 

3 

4 
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Routines Called: 
Output: 


In configuration 1 the minimum load value is sought within the limits of 
the low and high wave value. In configuration 2 the minimum is sought at lower 
wave values until the lower limit is hit or until a higher load value is found 
(which converts the configuration to type 1). In configuration 3 the minimum 
is sought at higher wave values until the upper limit is hit or until a higher load 
value is found (which converts the configuration to type 1). In the case of 
configuration 4, it is treated as a combination of configurations 2 and 3 and 
the lower of the two minimums from each configuration is found. 

There is no guarantee that the lowest minimum is found if more than one 
minimum can exist. (Multiple local minimums do occur in BUCLAP2 problems.) 

None 

(a) Calling sequence— none 

(b) COMMON— /NDIR/, NDIR 

/WAVE/, NOWAVE, WAVE, VLOAD 
/INFO/, SWAVE, KEND 

/WVSRCH/, NTRY, NFLIP, NHALF, BORDR1, BORDR2, 

BORDR3, WAVEP, WAVES, WAVET, VLOADP, 
VLOADS, VLOADT, HOLDL, HOLDW 
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Subroutine ZARK2 (N, MAX, EP1, EP2, FUN, I , ANSA, FANSA) 


Purpose: 
Called From: 
Input: 


Discussion: 


Routines Called 
Output: 


To find N zeros of an arbitrary complex valued function of a complex variable. 

AROOT 

Calling sequence 

N— number of zeros to be found 

MAX— maximum number of iterations to be executed in finding any one 
zero 

EP lb- convergence tolerances (Iteration will be terminated if |Xj - 
Ep 2 EP1. ( |Xj|) or |F (Xj) | < EP2. Xj and Xj.j are the last two appro- 
ximations to the zero.) 

FUN— complex function subprogram with complex parameter X which 
returns FUN = F (X) 

Muller’s method of successive approximation with quadratic complex polyno- 
mials. The Newtonian form of the approximating polynomial is used. Successive 
zeros are found by factoring out previously found zeros. That is, if F(Z) is the 
original function, and a zero Z Q is computed, then the next zero will be found 
for the function F(Z)/(Z - Z 0 ). In general, if Zj, i = 0, ..., n, are calculated zeros, 
the next zero will be sought for the function. 

F(Z)/ " (Z-Zj) 
i = 0 

The complex conjugate of a complex root is sought immediately. 

None 

Calling sequence 

ANS— N dimensional complex array returning the computed zeros 
FANS— N dimensional complex array returning the function values at the 
computed zeros. 

I— status flag 

= 0 run successful 

= 1 run unsuccessful— failed to converge in maximum iterations 
N— number of zeros computed 
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2.0 PROGRAM VALIDATION 


The program BUCLAP2 was validated by running various tests (data input sets unless noted 
otherwise) such that all major logic path and input options were tested. See reference 2 for data 
description. The validation testing was done in two ways: execution of actual data and inspection. 


2.1 INPUT OPTION FUNCTIONAL TESTS 


Card 

Field 

Option 

Test 

Cl 

1 

Data set title 

H-l 

C2 

4 

a root check print 

H-20 


5 

R matrix check print 

H-20 


6 

X^, Xf, stiffness check print 

H-20 


7 

Merged stiffness check print 

H-20 


8 

Load control check print 

H-20 


9 

Q print 

H-20 

C3 

1 

Plate upper bound only 

H-28 


2 

Date check only 

H-42 


3 

Number of buckling iterations 

H-10 

LC 

1 

Load case option 




= 1 

H-l 



2 

H-2 



3 

H-5 



4 

H-6 



5 

H-7 



6 

H-8 



7 

H-9 


2 

First load value— (non-blank) 

H-l 


3 

Second load value— (non-blank) 

H-6 


4 

Lower-bound load 

H-l 


5 

Upper-bound load 

H-l 
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Card 


Field 


Option 


Test 


W1 


W2 

W3 

BL 


1 

Half-wave input option 



= 1 

H-l 


2 

H-3 

2 

Loop control 



= 1 

H-10 


2 

H-l 


3 

H-l 1 


4 

H-5 


5 (wave search) 

H-3 

3 

Initial m or X values 

H-l 

4 

Final m or X values 

H-l 

5 

A m or A X values 

H-4 

1-15 

Wave number list 

H-l 1 

1 

Wave length table 

H-3 

1 

Preset boundary condition 



= SS 

H-l 


CL 

H-3 


FR 

H-29 

2 

= 1 

H-14 


2 

H-l 2 

3 

= 1 

H-14 


2 

H-12 

4 

= 1 

H-14 


2 

H-12 

5 

= 1 

H-30 


2 

H-14 

6 

= 0 

H-31 


k w 

H-12 

7 

= 0 

H-31 


k e 

H-12 

8 

= 0 

H-31 


k v 

H-12 

9 

= 0 

H-31 


k u 

H-12 
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Card 

Field 

Option 

Test 

BR 

1 

Preset boundary condition 




= SS 

H-l 



CL 

H-29 



FR 

H-3 


2 

= 1 

H-I4 



2 

H-l 2 


3 

= 1 

H-l 4 



2 

H-l 2 


4 

= 1 

H-l 4 



2 

H-l 2 


5 

= 1 

H-30 



2 

H-l 4 


6 

= 0 

H-31 



1c 

H-l 2 



w 



7 

= 0 

H-31 



k 0 

H-l 2 


8 

= 0 

H-31 



k v 

H-l 2 


9 

= 0 

H-31 



k w 

H-l 2 

TT 

1 

Table entries 

H-3 

AT 

1 

Table entries 

H-3 

MT 

1 

Table number 

H-3 


2 

Table entries 

H-3 


1 

Plate type 



= 1 

H-l 


2 

H-4 

2 

Number of laminas 

H-l 

3 

Input format 



= 0 

H-l 


1 

H-25 
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Card 


Field 


Option 


Test 


4 

Input option 



= 0 

H-l 


1 

H-3 

5 

Element length 

H-l 

6 

Element width 

H-l 

7 

Element radius 

H-4 




1 

Lamina thickness 

H-l 

2 

E 11 

H-l 

3 

E 22 

H-2 

4 

V 12 

H-l 

5 

G i2 

H-l 

6 

d 

H-3 2 
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Card 

Field 

Option 

Test 

P 3 

First set 

Thickness pointers 

H-3 


Second set 

Material pointers 

H-3 


Third set. 

Angle pointers 

H-3 

LA 

1 

Table entries 

H-83 

LM 

1 

Table entries 

H-83 

LT 

1 

Table entries 

H-83 

END 

1 

End of data set card 

H-l 
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2.2 DIAGNOSTIC FUNCTIONAL TESTS 


Error number 

Test 

50 

H-62 

51 

H-38 

52 

H-39 

53 

H-40 

54 

H-37 

60 

H-35 

61 

H-67 

62 

H-68 

70 

By inspection 

71 

By inspection 

75 

H-63 

80 

By inspection 

81 

By inspection 

82 

H-41 

83 

By inspection 

90 

By inspection 

91 

By inspection 

92 

By inspection 

95 

By inspection 

101 

By inspection 

102 

H-62 

105 

H-69 







2.3 PROGRAM FUNCTIONAL TESTS 


Program functional tests are tests of possible paths through each subroutine. They are 
listed in an alphabetical order of subroutines. 


Routine 

Configuration 

Test 

AROOT 

1 . Flat plate with B = 0 

H-l 


2. Flat plate with B^O 

H-2 


3. Check print of ot root information 

H-20 


4. a root which is real only 

H-74 


5. a root which is pure imaginary 

H-74 


6. a root which is arbitrary complex 

H-74 

ASTAR 

See description 

H-l 

CALVET 

See description 

H-l 

CDTM 

See description 

H-l 

CHMDET 

See description 

H-l 

DATIN 

1. Cl card present 

H-l 


2. Cl card not present - error 

H-33 


3. C2 card present 

H-l 


4. C2 card not present 

H-3 


5. C3 card present 

H-l 


6. C3 card not present 

H-3 


7. LDOPT = 1 

H-l 


8. LDOPT = 2 

H-2 


9. LDOPT = 3 

H-5 


10. LDOPT = 4 

H-6 


11. LDOPT = 5 

H-7 


12. LDOPT = 6 

H-8 


13. LDOPT = 7 

H-9 


14. IWTYPE = 1 

H-l 


15. IWTYPE = 2 

H-3 


16. W3 cards present 

H-3 


17. BL card with preset boundary conditions 

H-l 


18. BL card without preset boundary conditions 

H-14 
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Routine 


Configuration 


Test 



19. BR card with preset boundary conditions 

H-l 


20. BR card without preset boundary conditions 

H-l 4 


21 . Left side SS boundary condition 

H-l 


22. Left side CL boundary condition 

H-3 


23. Left side FR boundary condition 

H-29 


24. Right side SS boundary condition 

H-l 


25. Right side CL boundary condition 

H-29 


26. Right side FR boundary condition 

H-3 


27. Lamina thickness table 

H-3 


28. Material table 

H-3 


29. Angle table 

H-3 


30. Flat plate 

H-l 


31. Curved plate 

H-3 

DBLERT 

See description 

H-63 

\ 

DC 

1. LDOPT = 1 

H-l 


2. LDOPT = 2 

H-3 


3. LDOPT = 3 

H-5 


4. LDOPT = 4 

H-6 


5. LDOPT = 5 

H-7 


6. LDOPT = 6 

H-8 


7. LDOPT = 7 

H-9 


8. Check print of merge matrix 

H-20 


9. Matrix positive definite 

H-20 


10. Matrix not positive definite 

H-20 

DTC 

1 . Flat plate with B = 0 

H-l 


2. Flat plate with Bs^O 

H-2 

ERROR 

See description 

H-7 5 

GALUP 

1 . Flat plate 

H-l 


2. Curved plate 

H-3 


3. LDOPT = 1 

H-28 


4. LDOPT = 2 

H-3 


5. LDOPT = 3 

H-5 


6. LDOPT = 4 

H-6 
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Routine Configuration Test 

7. LDOPT = 5 H-7 

8. LDOPT = 6 H-8 

9. LDOPT = 7 H-9 

INFORM 1. Flat plate H-l 

2. Curved plate H-3 

3. LDOPT = 1 H-l 

4. LDOPT = 2 H-2 

5. LDOPT = 3 H-5 

6. LDOPT = 4 H-6 

7. LDOPT = 5 H-7 

8. LDOPT = 6 H-8 

9. LDOPT = 7 H-9 

10. IWTYPE = 1 H-l 

11. IWTYPE = 2 H-3 

12. Edge with SS boundary condition H-l 

13. Edge with CL boundary condition H-3 

14. Edge with FR boundary condition H-3 

15. Edge with input stiffness value H-2 5 

LOADCL 1 . Upper bound input H-l 

2. Upper bound calculated H-3 

3. IWTYPE =1 H-l 

4. IWTYPE = 2 H-3 

5. MSOPT = 1 H-10 

6. MSOPT = 2 H-l 

LOADCN 1 . First pass H-l 

2. Second pass H-l 

3. 0/1 crossing and determinant exponents < 8 H-7 6 

4. 0/1 crossing and determinant exponents >8 H-10 

MERGEC See description H-l 
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Routine 

NEXTC 

PAC 

PLATEC 

PLTDEF 

PLTCLC 


Configuration Test 

1 . Illegal character H-66 

2. Process new card H-3 

3. Process old card H-3 

See description Inspection 

1 .Flat plate with B = 0 H-l 

2. Flat plate with B # 0 H-2 

1 . Input by engineering constants by P2 cards H-l 

2. Input by engineering constants by MT cards H-3 

3. Input by matrices H-2 5 

4. Angle input on same card H-3 2 

5. Angle input separately H-3 

6. Flat plate B = 0 H-l 

7. Flat plate B * 0 H-2 

8. Curved plate H-3 


1 . Flat plate B - 0, 1,6 and 


2, 6 zero 

AN1 1 = 0 

AN12 = 0 

AN22=£0 

H-17 

2. 


= 0 

#0 

= 0 

H-10 

3. 


# 0 

= 0 

= 0 

H-l 6 

4. 


# 0 

#0 

= 0 

H-l 1 

5. 


= 0 

¥=0 

=£0 

H-l 2 

6. 


#0 

= 0 

¥=0 

H-l 

7. ’ 

f 

#0 

#0 

=£0 

H-6 

8. Flat plate B # 0 

= 0 

= 0 

#0 

H-77 

9. 


= 0 

=£0 

= 0 

H-78 

10. 


#0 

= 0 

= 0 

H-79 

11. 


#0 

=# 0 

= 0 

H-80 

12. 


= 0 

=£0 

¥=0 

H-81 

13. 


#0 

= 0 

# 0 

H-2 

14. 1 


=£ 0 

#0 

#0 

H-82 

15. Curved 

jlate B = 0 

= 0 

= 0 

# 0 

H-21 

16. 


= 0 

¥=0 

= 0 

H-13 

17, 


* 0 

= 0 

= 0 

H-20 

18. f 

=£0 

# 0 

= 0 

H-l 5 
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Routine 


Configuration Test 


19. Curved plate B = 0 

AN1 1 =0 

AN 12 #0 

AN22 # 0 

H-23 

20. | 


=£0 

= 0 

=£0 

H-22 

21, 1 

f 

¥= 0 

=£0 

¥= 0 

H-24 

22. Curved plate B ¥= 0 

= 0 

= 0 

# 0 

H-l 8 

23. 


= 0 

#0 

= 0 

H-l 4 

24. 


0 

= 0 

= 0 

H-l 9 

25. 


#0 

^0 

= 0 

H-25 

26. 


= 0 

0 

#0 

H-26 

27. 


¥= 0 

= 0 

¥= 0 

H-4 

28, j 


^ 0 

=£ 0 

0 

H-27 


29. Curved plate B 0, 1,6 and 


2, 6 terms nonzero H-5 

30. Flat plate B = 0, 1,6 and 

2, 6 terms nonzero H-3 

PRERD 1 . No end card H-34 

2. Cl card present H-l 

3. C2 card present H-l 

4. C3 card present H-l 

5. LC card present H-l 

6. W.l card present H-l 

7. W2 card present H-l 1 

8. W3 card present H-3 

9. BL card present H-l 

10. BR card present H-l 

1 1 . TT card present H-3 

12. MT card present H-3 

13. AT card present H-3 

14. PI card present H-l 

15. P2 card present H-l 

16. P3 card present H-3 

17. PA card present H-25 

18. PB card present H-25 

1 9. PD card present H-25 

20. END card present H-l 

See description H-20 


PRINTC 
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Routine 

Configuration 

Test 

RDTBLE 

1. Two or more cards used 

H-36 


2. E format data 

H-3 


3. E format data with 14 significant digits 

H-36 


4. F format data 

H-3 


5. E format data with 14 significant digits 

6. Error exits 

H-36 


(a) Missing slash 

H-3 7 


(b) Illegal data (bad character or too many characters) 

H-3 8 


(c) Table size limit exceeded 

H-3 9 


(d) No data found 

H-65 

REDUC3 

1 . Successful operation 

H-3 


2. Unsuccessful operation 

H-41 

RGENER 

Same as for PLTCLC 

Same as PLTCLC 

SELECT 

1 . Wave search option 

H-3 


2. Non-wave search option 

H-l 


3. Upper bound only 

H-28 


4. Upper bound input 

H-l 


5. Full solution 

H-l 

SOLVEC 

See description 

H-l 

STOREC 

See description 

H-l 

STRAIN 

1, LDOPT = 1 

H-l 


2, LDOPT =2 

H-3 


3. LDOPT = 3 

H-5 


4. LDOPT = 4 

H-6 


5. LDOPT = 5 

H-7 


6. LDOPT = 6 

H-8 


7. LDOPT = 7 

H-9 

STRMOV 

See description 

Inspection 

S0352A 

1 . Full solution 

H-l 


2. Upper Bound only 

H-28 
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Test 

H-l 

H-2 

H-3 

H-4 

H-5 

H-6 

H-7 


2.4 FUNCTIONAL TEST DESCRIPTIONS 


Characteristics 


1 . Flat plate, matrix B = 0 

2. Cl title card used 

3. Load option = 1 , input N y 

4. BL = BR = SS 

5. Input upper bound 

6. Wave number option 

7. MSOPT = 2 

8. Input engineering constants 

9. END card included 

1 . Flat plate, matrix B # 0 

2. Load option = 2, input N x 

1 . Curved plate with matrix B = 0 

2. Load option = 2 

3. BL = clamped 

4. BR = free 

5. Upper bound defaulted 

6. TT, MT and AT tables used 

7. Minimum wave search option used— MSOPT = 5 

8. Half wave length option 

9. C2 and C3 cards not present 

1. Curved plate, matrix B # 0; Bjg and B 26 = 0 

2 . Load option = 2, input N x 

1. Curved plate, matrix B 0, Bjg and 625 ^ 0 

2. Load option = 3 

3. MSOPT = 4 

1 . Flat plate, matrix B = 0 

2. Load option = 4, input N x and N y /N xy 

1 . Flat plate, matrix B = 0 

2. Load option = 5, input N y and N x /N xy 
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Routine 

Configuration 

Test 

TBPOINT 

1 . Repeat option left of slash 

H-3 


2. Repeat option right of slash 

H-3 

UNPAC 

See description 

Inspection 

VIPDR 

See description 

H-l 

WAVFND 

1 , Load value for initial middle wave value is lowest of the three 

H-70 


2. Load values corresponding to initial three wave values are 
increasing 

H-71 


3. Load values corresponding to initial three wave values are 
decreasing 

H-72 


4. Load value for initial middle wave value is highest of the three 

H-73 

ZARK2 

See description 

H-l 
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Test 


Characteristics 


H-8 

H-9 

H-10 
H-l 1 

H-12 

H-l 3 
H-l 4 
H-l 5 
H-l 6 
H-l 7 
H-l 8 


1 . Flat plate, matrix B = 0 

2. Load option = 6, input N x /N y 

1 . Flat plate, matrix B = 0 

2. Load option = 7, input N x /N y and N y /N xy 

1 . Flat plate, matrix B = 0 

2. Load option = 3 

3. MSOPT = 1 for wave number loop control 

4. Buckling iterations not defaulted 

1 . Flat plate, matrix B - 0 

2. Load option = 3, input N x 

3. MSOPT = 3 

4. W2 card used 


1 . Flat plate, matrix B = 0 

2. Load option = 2, input N x 

3. Spring constants used for boundary 


1 . Curved plate, matrix B = 0 

2. Load option = 3 


1 . Curved plate, matrix B # 0 
2 . Load case = 3 

1 , Curved plate, matrix B = 0 

2. Load case = 3, input N x 

1 . Flat square plate, matrix B = 0 

2. Load option - 1 

1 . Flat plate, matrix B = 0 

2. Load option - 2 

1 . Curved plate, matrix B # 0 

2 . Load option = 2 



Test 

H-19 

H-20 

H-21 

H-22 

H-23 

H-24 

H-25 

H-26 

H-27 

H-28 

H-29 

H-30 


Characteristics 

1 . Curved plate, matrix B =£ 0 

2. Load option = 1 

1 . Curved plate, matrix B = 0 

2. Load option = 1 

3. Intermediate check print 

1 . Curved plate, matrix B = 0 

2. Load option = 2 

1 . Curved plate, matrix B = 0 

2. Load option = 1 , input N y 

1 . Curved plate, matrix B = 0 

2. Load option = 3, input N y 

1. Curved plate, matrix B = 0 

2. Load option = 7, input N x /N y and N y /N xy 

1 . Curved plate, matrix B ¥= 0 

2. Load option = 5, input N x /N xy 

3. Input A, B and D matrices 

1 . Curved plate, matrix B =£ 0 

2. Load option = 4, input N y /N xy 

1 . Curved plate, matrix B ¥= 0 

2. Load option = 7, input N x /N y and N y /N xy 

1 . Same plate as H-l 

2. Calculate upper bound only 

1 . BL = free 

2 . BR = clamped 

1 . Boundaries restrained 
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Test 

H-31 

H-32 

H-33 

H-34 

H-35 

H-36 

H-37 

H-38 

H-39 

H-40 

H-41 

H-42 


H-43 


Characteristics 

1 . Free boundaries 

1 . Fiber angle input by P2 card 

1. Cl card not present 

1 . No END card 

1. Illegal character in P3 card 

1 . Two cards used for W3 card 

L Slash missing in W3 card 

1 . Illegal character in W.3 card 

1. W3 card table size limit exceeded 

1. No data found on W3 card 

1. Size load greater than critical load, plate 
already buckled 

1 . Stack of H3 data sets 

2. LC card missing in first set 

3. W1 card missing in 2nd set 

4. W3 card missing in 3rd set 

5. BL card missing in 4th set 

6. BR card missing in 5th set 

7. TT card missing in 6th set 

8. MT card missing in 7th set 

9. AT card missing in 8th set 

10. PI card missing in 9th set 

1 1 . P3 card missing in 10th set 

1 . H-2 except JPC(3) = -10 on C3 card 
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Characteristics 


Test 

H-44 

H-45 

H-46 

H-47 

H-48 

H-49 

H-50 

H-51 

H-52 

H-53 

H-54 

H-55 

H-56 

H-57 

H-58 


1 . H-2 except load option = 8 

1. H-3 except 41 material tables input 

1 . H-3 except plate type = 0 

1 . H-3 except material property option = 2 

1 . H-3 except material properties input format 
option = 2 

1. Problem with 26 layers 

1 . H-3 except wave option = 1 used with initial 
half-wave length value higher than second and a 
negative delta supplied 

2, Half-wavelength list exceeds 30 values 

1 . H-3 except BL = BL 

1. H-3 except wave-value definition on the W1 card 
is left blank 

1 . H-3 except wave value definition is set equal to 1 

1. H-3 except for the following items: 

2. MT table has 3 values 

3. Radius is set negative 

4. Length value is given 

1 . Insufficient P.2 cards 
1 . Invalid table number 
1 . An input E 22 value is zero 
1 . An input E j j value is zero 
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Test 


H-59 


H-60 

H-61 

H-62 

H-63 

H-64 

H-65 

H-66 

H-67 

H-68 

H-69 

H-70 

H-71 

H-72 


Characteristics 

1 . Data set containing a card which is an invalid 
type 

1 . Data set with P2 and P3 cards 

1 . Data set with a missing P3 card 

1. An end-of record card following an MT card 

1 . Double root occurs 

1 . An MT card starts list with a comma 

1 . W3 card with a slash only 

1 . Illegal character on a P3 card 

1 . P3 card starts with 2 commas 

1 . Unequal number of commas in P3 sublists 

1 . Upper bound input as too low 

2. Critical wave search option chosen 

1 . Critical wave search option is input 

2, Wave-load configuration 1 (see subroutine 
description) 

1 . Critical wave search option is input 

2. Wave-load configuration 2 (see subroutine 
description) 

1 . Critical wave search option is input 

2. Wave-load configuration 3 (see subroutine 
description) 
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Test 

H-73 

H-74 

H-75 

H-76 

H-77 

H-78 

H-79 

H-80 

H-81 

H-82 

H-83 


Characteristics 

1 . Critical wave search option is input 

2. Wave-load configuration 4 (see subroutine 
description) 


1. a roots of each type— pure real, pure imaginary, 
and complex— occur. (This is not an input data set; 
subroutine ZARK2 altered for the test.) 

1 . Error exit through subroutine ERROR occurs 

1 . Test of LOADCN with determinant exponents < 8 

1. B ¥= 0 

2. Nil = 0 and N1 2 = 0 

3. N22 - variable 

1. B¥=0 

2. N1 1 = 0 and N22 = 0 

3. N1 2 = variable 

1. B ¥= 0 

2. N12 = 0 and N22 = 0 

3. N1 1 = variable 

1. B ¥=0 

2. N1 2= 1000,0 and N22 = 0 

3. Nil = variable 

1. B ¥= 0 

2. N1 1 = 0 and N12 = 1000.0 

3. N22 = variable 

1. B ¥= 0 

2. N1 1 = 1000.0 and N12 = 60.0 

3. N22 = variable 

LT, LA, LM cards are used in data input. 
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ENGINEERING VARIABLE GLOSSARY 


This shows the correspondence between the engineering variables in the analysis document, 


reference 1, and the program variables. 
Engineering 

Program 

variable 

Description 

variable 

a 

length of the plate 

PLONG 

A ij 

extensional stiffnesses (i, j = 1 , 2, 6), equation (A-6) 

A 

A* 

matrix for strain calculation, equation (A-56) 

AST 

b 

developed width of the plate 

WIDE 

B ij 

stiffnesses (i, j = 1, 2, 6) associated with bending- 
stretching coupling, equation (A-6) 

B 

D ij 

bending stiffnesses (i, j = 1 , 2, 6), equation (A-7) 

D 

E 11’ E 22 

Young’s moduli of ortho tropic lamina 

El, E2 

G 12 

shear modulus of orthotropic lamina 

G12 

h k 

distance to lamina from the reference plane 

H(K) 

k^y* k^ 

spring constants 

DSTIFF 

*1. *2 

direction cosines 

X, Y 

E 1 y L 2j 

displacement ratio coefficients, equations (A-38) 
and (A-39) 

L1,L2 

ill 

moment resultant, equation (A-15) and (A- 19) 

CZ 

n 

number of laminas 

LAYERS 

N N N 

stress resultants, equation (A-6) 

STR 
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Engineering 


Program 

variable 

Description 

variable 

N X , Ny, N x y 

applied inplane loads 

AN1 1, AN12, AN22 


buckling displacement parameter 

P 

Q ij 

orthotropic material constants (i, j = 1,2, 6) with 

Q1 1, Q12, Q22, 


respect to material axes, equation (A-5) 

Q33 


orthotropic material constants (i, j = 1 , 2, 6) with 
respect to plate axes, equation (A-3) 

Q 

A 



Q 

effective transverse shear, equations (A-14) and (A-18) 

cw 

R 

reference plane radius of curved plate 

RADIUS 

R 1 1» R 12» • • • R 33 

elements of coefficient matrix R, equation (A-37) 

R 

s 

stiffness matrix for the plate, equation (A-53) 

ST 

S 

modified merged stiffness matrix, equation (A-54) 

DBMA 

l k 

thickness of k^ layer of the laminate 

TH(K) 

t 

laminate thickness 

TOT 

A 



T 

effective inplane shear, equations (A- 17) and (A-20) 

CU 

x d ,x f 

displacement and force matrices, respectively, 
equations (A-45) and (A-5 1 ) 

XD, XF 


buckling displacement parameters, equations (A-34) 
through (A-36) 

ALPHA, BETA 

Txy 

shear strain 

STR(3) 

e x , e y 

normal strains 

STR(l), STR(2) 
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Engineering 

variable 


Description 


Program 

variable 


X 

v \V v 2\ 

^k 

Subscripts 


half-wavelength of buckling in the x-direetion SWAVE 

Poisson’s ratio RNUA, RNUB 

angle defining orthotropy directions of k^ 1 lamina, ANGLE 
with respect to plate axes 


j index corresponding to characteristic roots 

k layer index 

A subscript preceded by a comma indicates partial differentiation with respect to the subscript. 
Superscripts 


T 


o 


+ 


matrix transpose 

quantities in the reference plane of the plate 
quantities along the side y = of the plate 
quantities along the side y = ~ of the plate 
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